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Buchert et al have now published a revised version [T] of their original paper [2] 
criticizing our work. In the revised version, the tone of their criticisms has been improved 
and all assertions that we made mathematical errors in our work have been removed. 
However, the essential content of their paper remains unchanged. Therefore, we shall 
hereby revise our rebuttal by simply adding some remarks to our original rebuttal to 
note changes and address new points. In addition, elsewhere [3] we will give a simple, 
heuristic discussion of why backreaction is negligible in cosmology. Except for the 
remarks identihed below as ^^[NOTE ADDED: . the remainder of this paper is 
unchanged from our original posting. 

In a recent paper, Buchert et al [2] have criticized nearly all aspects of our analysis 
of backreaction produced by small scale inhomogeneities in cosmology as presented in 
[7j. Most of these criticisms concern points that are irrelevant to our main arguments 
and/or are based upon misinterpretations of our work. Nevertheless, in order to avoid 
the appearance that we are ignoring or evading potentially valid issues, we will give 
a brief response in the Appendix to all of the scientihc/mathematical points raised 
in [2]. In short, we stand hrmly behind all of the assertions, arguments, and conclusions 
presented in all of our previous papers iHZl. 

Nevertheless, we believe that the discussion of [2] has the potential to create some 
confusion about the nature of our analysis of backreaction and the conclusions that 
can be drawn from it. Therefore, we feel that it would be appropriate for us to take 
this opportunity to explain, in an informal way, what we mean by “backreaction” and 
how the main tools of our analysis £t in with prior efforts, particularly the work of 
Isaacson Eli- 

The context for the phenomenon of backreaction by small scale inhomogeneities 
concerns a situation where the actual spacetime metric Qab has large curvature 
fluctuations on small scales but, nevertheless, Qab can be well approximated by a metric 

that does not have large curvature fluctuations. Although our analysis is valid in 
a much more general context, the main situation we have in mind is where gab is the 
actual metric of the universe and g^J is a metric with FLRW symmetry. For simplicity, 
we will restrict our discussion below to the case where g^^fj has FLRW symmetry. 

The issue at hand is whether the small scale inhomogeneities of gab can contribute 
nontrivially to the dynamics of g^^fj. A priori, this is possible even though 'jab = Qab — g^ab 
is assumed to be small: Einstein’s equation for gab contains derivatives of 7 ^ 6 , which 
need not be small even when 7^6 is small. Consequently, the Einstein tensor. Gab, of gab 
need not be close to the Einstein tensor, Gjll\ of gj^^K Thus, although gab is assumed 
to be an exact solution of Einstein’s equation with some stress-energy source Tab, h is 
possible that may not be close to a solution to Einstein’s equation with a suitably 
averaged stress-energy source If this occurs, we say that there is a substantial 

backreaction effect of the small scale inhomogeneities on the effective dynamics of g'^^^. 

There are two cases where our analysis does not apply (and was never intended or 
claimed to apply). The hrst case is where the actual metric gab is not close to any FLRW 
metric g^^f^. A good example of this relevant to cosmology is obtained by taking gab to 
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be an LTB model [IQHII] with significant voids/overdensities on scales comparable to 
the Hubble radius. A simpler (but much less physically motivated) example discussed 
in na is the case where the universe consists of two disconnected FLRW parts at 
different stages of expansion/contraction. In such cases, since gab differs substantially 
from any single FLRW metric, there is little point in even attempting to hnd a “best £t” 
FLRW metric by averaging or by any other scheme. If one insists on hnding a “best 
£t” g^J, the dynamics of g^^J would be expected to differ signihcantly from standard 
FLRW dynamics. However, we would not describe this as “backreaction.” 

The second case to which our analysis does not apply is where the actual metric 
gab is close to an FLRW metric g^^J, but one applies an averaging or other scheme to 
gab to construct an “effective” FLRW metric intended to reproduce averages of 
observables calculated in the actual metric. In general, this approach will give rise to a 
metric that differs signihcantly from g^J. This will typically happen if one tries to 
match the geodesics of gab in a naive way to corresponding geodesics in an FLRW model, 
since—although gab is close to g'fj ^—the geodesics of gab are not close to the geodesics 
of g^^fj. We provided an illustration of this in our “ball bearing” example in [7], where 
we constructed a smooth metric Qab on a 2-sphere that was everywhere close to the 
round sphere metric but the geodesics and curvature of Qab are not close to those 
of As we explained in |7j, observers living in “Sphereland” and doing experiments 
with taut ropes (geodesics) would hnd extremely large deviations from any round sphere 
model on small scales and also would hnd some signihcant deviations on large scales. If 
they used their geodesic and curvature data to determine a “best ht” round metric q}l^\ 
they could easily make a substantial error, depending upon exactly what “averaging 
scheme” they used. Returning to the cosmological case, we noted in section 3 of [7] 
that the Buchert approach uses timelike geodesics in an essential way and the Clarkson 
and Umeh approach uses null geodesics in an essential way, so a priori both approaches 
may produce an FLRW metric that dihers signihcantly from g^^^. The dynamics of 

may then diher signihcantly from that of If this occurs, we shall refer to the 
ehect as pseudo-backreaction, since it has the appearance of being a dynamical ehect 
produced by small scale inhomogeneities but is actually merely an artifact produced 
by a poor choice of representative FLRW metric Our analysis does not apply to 
pseudo-backreaction. 

If, as we claim, naive averaging schemes making use of geodesics and curvature are 
prone to producing errors in the choice of representative FLRW model g'^J, how should 
one go about determining As we discussed in section 4.3 of [7j, this should be done 
in the manner that is normally already done in practice by cosmologists (see, e.g., US]): 
Consider FLRW metrics (with appropriate free parameters) and their perturbations 
(with appropriate free parameters). Calculate observable quantities within the models 
and hnd the model parameters that provide the best ht to all of the data. If one does 
not hnd an acceptable ht, this shows that something is wrong or incomplete in the 
theoretical model. If one does hnd an acceptable ht, this determines the parameters 
of the model, and one thereby obtains a representative FLRW metric g^^fj. Of course. 
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this does not “prove” that the theoretical model is correct or that is close to the 
actual metric of the universe, Qab- As in all other areas of science, even if one has a 
simple model that provides an extremely good £t to a wide variety of disparate data, 
one always should be open to the possibility that there are alternative explanations. 

None of what has been said above addresses the issue of whether, in our actual 
universe, there might be signihcant backreaction effects on the large scale dynamics 
produced by small scale inhomogeneities. As discussed above, this is possible a priori 
because even if ^ah = Qab — 9^ab small, its contribution to Einstein’s equation need not 
be small. How can one compute these backreaction effects? 

The issue of backreaction was addressed in 1964 by Brill and Hartle ca in the 
context of trying to hud “geon solutions” to Einstein’s equation, consisting of a ball 
of gravitational radiation that is held together by its self-gravitation. In this case, 
the small scale inhomogeneities consist of gravitational radiation, whose backreaction 
substantially alters the background metric. The Brill-Hartle approach was then 
signihcantly generalized by Isaacson mu- The Isaacson work is nicely summarized 
in subsections 35.13-15 of Misner, Thorne, and Wheeler na. 

However, although the approximations made in Isaacson’s work are well motivated 
physically, some serious difficulties arise if one tries to give a mathematically precise 
justihcation of them. The most serious difhculties involve the justihcation of equations 
satished by ■jab (denoted as in [IE]), such as eq. (35.59a) of [TH]. The nature of these 
difficulties were elucidated at the end of section HI of our hrst paper |1]. 

In the late 1980’s, one of us (R.M.W.) was very troubled by these difficulties 
and felt that the best way to resolve them would be to reformulate the Isaacson 
approximation scheme in terms of one-parameter families of metrics gab{)^), similar 
to what is done to justify ordinary perturbation theory (see section 7.5 of PS]). A 
reformulation of the Isaacson scheme in terms of one-parameter families would allow 
one to say with much more clarity what approximations are being made and would 
enable one to rigorously derive the equations satished by the various quantities at 
each order of approximation. However, to describe the Isaacson scheme in this way, 
one would need to consider one parameter families where 7 a 6 (A) —)■ 0 as A —>■ 0 
but spacetime derivatives of 7 a 6 (A) do not go to zero, in order that that they may 
continue to make nontrivial contributions to the second order Einstein tensor in the 
limit, thereby providing nontrivial backreaction. It is therefore far from obvious how to 
describe the limiting behavior as A —)■ 0 in a mathematically precise way. G. Burnett 
was a graduate student in the Chicago relativity group and R.M.W. posed to Burnett 
the problem of reformulating the Isaacson approximation scheme using one-parameter 
families. Burnett solved this problem brilliantly in [20], using the notion of weak limits 
to give a mathematically precise characterization of the one-parameter families gab{)^)- A 
discussion of the relationship between Burnett’s reformulation and the original Isaacson 
scheme can be found in Burnett’s paper and our previous papers (particularly sections 
I and H of jl]). 

Our contribution to the analysis of backreaction was to recognize that the 
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methods used in the Isaacson approach—as reformulated by Burnett—could be directly 
imported to describe backreaction in cosmology. The small scale inhomogeneities of 
primary interest now are not gravitational waves oscillating in space and time but 
are the perturbations associated with density inhomogeneities varying mainly in space. 
Nevertheless, the scheme for calculating the “Isaacson average” of the second order 
Einstein tensor of these perturbations does not depend on any wavelike character 
of the perturbation and works just as well for treating the backreaction effects of 
Newtonian-like cosmological perturbations as it does for treating the backreaction effects 
of gravitational waves. 

Our general analysis presented in |1] allowed for the presence of gravitational 
radiation. It thereby allowed for backreaction produced by gravitational radiation as well 
as by matter inhomogeneities. Our main results are contained in two theorems. These 
theorems state that if the matter stress-energy tensor satisfies the weak energy condition, 
then the effective stress-energy tensor describing the leading order backreaction effects 
of the small scale inhomogeneities must (i) be traceless and (ii) satisfy the weak energy 
condition. In essence, our theorems say that significant backreaction effects can be 
produced only by gravitational radiation and not by matter inhomogeneities. 

At sub-leading order, matter inhomogeneities do produce nontrivial backreaction 
effects, as we analyzed in [^. In particular, in Appendix B of [5] we computed the 
backreaction effects of matter inhomogeneities on the expansion rate of the universe, 
under the assumption that these inhomogeneities are Newtonian-like in nature. A 
lengthy calculation revealed that backreaction effectively modifies the matter stress- 
energy by adding in the effects of kinetic motion and the Newtonian potential energy 
and stresses. In particular, it “renormalizes” the proper mass density to an “ADM mass 
density.” 

Our results are fully in accord with the “back of the envelope” estimates of 
backreaction given in na. Onr results on the sub-leading backreaction effects of 
Newtonian-like matter perturbations are fully in accord with the analysis of [21]. 
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Appendix A. Response to Specific Criticisms 

In this Appendix, we will provide brief responses to the specific criticisms of [2]. 
However, we will not address the following two categories of criticism: (1) Criticism 
or discussion of any results not presented in |1H7|; (2) criticisms that are of no scientific 
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importance. The second category includes the subcategories of (a) criticisms that are 
unr elated to the main issues at hand and (b) numerous insinuations and innuendo^ 
that suggest that we were sloppy in our assertions and arguments and/or overlooked 
various points. [NOTE ADDED: The tone has been improved in m-j Our silence 
on these categories of criticism should not be interpreted as indicating agreement. The 
section numbers below correspond to the numbering of [2]. 

Section 2: This section criticizes our ball bearing example (mentioned in the body of 
the paper above) where we constructed a smooth metric qab on a 2-sphere that was 
everywhere close to the round sphere metric but the geodesics and curvature of qab 
are not close to those of The purpose of this example was simply to illustrate how 
two metrics could be close without their geodesics and curvature being close. Therefore, 
logically, it does not make sense to criticize the example unless we made an error in 
our claim that the metrics are close or our claim that their geodesics and curvature are 
not close. Neither claim appears to be challenged in section 2. Thus, as a matter of 
principle, we do not understand what is being criticized. 

Section 3: This section points out that the Buchert formalism violates the conclusions 
of our theorems on the effective stress-energy tensor. The reasons why our results do 
not apply to the Buchert formalism were already explained in the body of the paper 
above. The remainder of this section mainly reviews our approach and raises numerous 
mathematical questions concerning the validity of our work, referring to Appendix B for 
details. We will give our response to these points under “Appendix B” below, except for 
the following comment: In contrast to what is said in subsection 3.7, backreaction in our 
formalism does not suddenly “turn on” at A = 0; the one-parameter family of metrics 
gabW converges uniformly on compact sets to so the dynamics of g^^J accurately 
describes the (large scale) dynamics of gabW for sufficiently small A. Indeed, for A > 0, 
backreaction terms are present, and are described to leading order by the second order 
Einstein tensor; for A = 0, they are accounted for by the effective stress-energy tensor 
t^^l\ which is the weak limit as A —)■ 0 of the second order Einstein tensor. 

Section 4 '■ Subsection 4.1 makes the point that whenever one has a one-parameter family 
of metrics 5 'ab(A), one could apply a A-dependent diffeomorphism to generate a new 
one-parameter family of metrics 'ip*{X)gabW- If 9abW approaches a limit as A —)■ 0, then 
a '4’{X) may be chosen so that ^p*{X)gab{^) approaches a different limitj§]or no limit at 
all. In this sense, the limit of any one-parameter family of metrics (or any other tensor 
fields) is always “coordinate dependent.” The situation with regard to our approach is 
no different from ordinary perturbation theory (see section 7.5 of |T9]) where one could 

f As a small but typical example, footnote 5 of [2] states that “Green and Wald obviously use a 
fixed A-independent Riemannian metric for computing norms...” The word “obviously” suggests that 
we omitted to say this. In fact, we stated that we were using a fixed, A-independent Riemannian 
metric for computing norms very clearly and explicitly in our papers. [NOTE ADDED: The word 
“obviously” has been deleted from footnote 5 in m-] 

§ The fact that different limits can be obtained in this manner for spacetimes representing a body that 
“shrinks to zero size” plays a central role in the Gralla-Wald derivation of self-force [221 US] ■ 
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start with a one-parameter family gab{)^,x) that is jointly smooth in A and x, and then 
apply a diffeomorphism 'ip{X) that is not smooth in A at A = 0 (as in the examples 
given in snbsection 4.1) so that 'il)*{X)gab{X) approaches a different limit (or no limit) 
as A —)■ 0. Obviously, the fact that one can perform such a construction is completely 
irrelevant to the validity of ordinary perturbation theory, and it is similarly irrelevant 
to the validity of our results. 

A much less trivial issue—which does not appear to have been considered in [2]— 
concerns the case where 'ip*{X) goes to the identity as A —)■ 0 (so that w-limA^-o = 
dlib) derivatives of 'i/j*{X) do not go to zero. Such transformations are of no 
relevance if 'ip*{X)gab{X) fails to satisfy our conditions (i)-(iv), just as non-jointly- 
smooth diffeomorphisms are of no relevance to ordinary perturbation theory. However, 
if 'ip*{X)gab{X) continues to satisfy our conditions (i)-(iv), it gives rise to a gauge 
dependence of [labcdef = w-limA-s-o Va7bcVd7e/ analogous to the gauge dependence 
of quantities appearing in ordinary perturbation theory. Nevertheless, as shown by 
Burnett m, the effective stress-energy tab describing backreaction is gauge invariant. 
This corresponds to the statement that the Isaacson average of the second order Einstein 
tensor becomes gauge invariant in the short-wavelength limit. Thus, our backreaction 
results are fully “coordinate invariant” in the desired sense. 

In subsection 4.2, our example given in [Hj of a family of Gowdy metrics exhibiting 
nontrivial backreaction is criticized on the grounds that it is “ultra-local” and also 
on the grounds that we could have chosen a different familj|j]] that would have given 
no backreaction. Our purpose in giving the Gowdy family was to provide an explicit 
example of a one-parameter family of metrics gab{X) that satishes our assumptions (i)- 
(iv) with nonvanishing backreaction, tab fy 0. Logically, it can be invalidated only by 
showing that our choice of gab{X) violates at least one of our assumptions (i)-(iv) or 
that tab = 0. This is not claimed/shown in subsection 4.2. 

Section 5\ This section is largely a defense against our previous criticisms of the Buchert 
approach. We stand by our previous criticisms. Subsection 5.4 discusses the issue of 
how to obtain the “smoothed” metric g^^^ from gab- As already discussed above in the 
body of the paper, it is our view that using geodesics and/or curvature of gab to directly 
construct g^J is a bad way to proceed, since it is prone to making signihcant errors 
because the geodesics and curvature of the actual metric gab differ signihcantly from 
gl^J. Rather, one should determine g^^J by building a model for g^J and 'jab and htting 
all available data to the model. 

Section 6: It is difficult for us to determine exactly what is being criticized in this section, 
but we take this opportunity to make two clarifying remarks related to the discussion: 
(1) Our dictionary given in [5] that translates a Newtonian model into a general 
relativistic spacetime will yield a spacetime metric that solves Einstein’s equation to 

II Note that our choice of one-parameter family with nonvanishing backreaction gives rise to a that 
provides a much better approximation to gabiXo) than the alternative families suggested in subsection 

4.2. 
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very high accuracy. The distinction between having a quantity nearly solve an equation 
at all times versus having this quantity provide a good approximation to a solution at 
all times was explained in the second-to-last paragraph of the introduction of [5]. (2) If 
gab is close to , then by flux conservation, the average apparent luminosity of sources 
(including multiple images) in the metric gab must closely match that of g^^f^. This is a 
simple fact, not an “argument against backreaction.” [NOTE ADDED: In some 
text has been added to try to rebut this simple fact. Here the authors confuse average 
apparent luminosity (which closely matches that of the background FLRW metric, by 
conservation of flux) with average luminosity distance (which is a nonlinear function of 
apparent luminosity and thus can differ significantly from an FLRW model). The fact 
that nonlinear functions of average apparent luminosity can be significantly affected by 
inhomogeneities was explained clearly in section f.3 ofE-l 

Appendix B: Appendix B contains the core of the arguments against our results. It 
is referred to many times in the body of |2] to justify their criticisms. Words such 
as “error” appear frequently in Appendix B when describing of our work. [NOTE 
ADDED: In these words have been removed.] 

The authors of [2] appear to assume that we wished to consider metrics gabW of 
low regularity in spacetime. This is not the case. Indeed, until we read [2j, it never 
occurred to us to consider spacetime metrics that are not C°° in spacetime at each fixed 
A, and there is no reason to consider metrics that are not smooth in spacetime. [NOTE 
ADDED: In m, the authors assert that “one cannot safely carry out computations by 
declaring that second derivatives of the perturbing tensor 7 may be unboundedly large, 
and then treat them as ordinary smooth functions. ” This suggests that some of the 
misunderstandings of the authors may originate from a confusion between unboundedness 
of second spacetime derivatives in the limit A —)■ 0 (which is allowed in our formalism) 
with unboundedness of these derivatives in spacetime at fixed A (which we had never 
contemplated allowing).] For smooth spacetime metrics, all of our manipulations are 
trivially justified. For example, in their eq. (B.9) one can simply integrate by parts 
a second time to take all the derivatives off of 7 ^ 6 , and then use our assumption (ii) 
(eq. (7) of |2]) to see that the limit vanishes. 

Nevertheless, the discussion of Appendix B raises an issue of some mathematical 
interest: If one did wish to consider metrics of low spacetime regularity, precisely what 
regularity would be needed to make our analysis valid? Indeed, since Valbc need not 
have a pointwise limit as A —)■ 0 and V flUblcd may become unboundedly large as A —0, 
it would be surprising if our results actually required smoothness of 7 ab(A) for all fixed 
A > 0. 

To analyze the regularity needed, we note that our condition (ii) (namely, 
| 7 afe(A,a;)| < ACi(a;)) makes sense only if 7^6 is a tensor field that is dehned (almost) 
everywhere. Condition (ii) then immediately further implies that 7afe is locally 
and thus defines a distribution. In particular, its (weak) second derivative^ are 

% [NOTE ADDED: What we are calling the “weak derivative” here corresponds to what the authors 
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automatically well defined as a distribution via 


= lcd[VbVar^^ = / laN.VJ 


cabcd 


(A.l) 


' M 


for any test (C'“) tensor field Similarly, condition (iii) (namely, | Va 7 fec(A, a;)| < 

C 2 {x)) makes sense only if 7 ^^ is differentiable (almost) everywhere. Condition (iii) 
then immediately further implies that Va^bc is locally L^. It follows immediately from 
the above together with the smoothness that the Einstein tensor of gab{^,x) = 
g'fb ix)+'yab{^, x) is well defined as a distribution on spacetime for all (sufficiently small) 
A. Thus, for all A > 0, it makes sense to require that Einstein’s equation holds weakly. 
Our results then follow by imposing Einstein’s equation weakly at A > 0 and taking the 
weak limit as A —)■ 0. The linear terms in the second derivatives of —the main focus 
of the critical remarks of Appendix B—are easily seen to make no contribution]^ in this 
limit on account of (A.l) and condition (ii). 


Thus, our results continue to hold under the weakest regularity conditions under 
which our assumptions make any sense. This, of course, is a side-point since we were, 
in any case, considering only metrics that are smooth in spacetime at each hxed A— 
in which case all of our manipulations are much more trivially justified. We conclude 
that the criticisms of Appendix B are completely without merit. Similarly, all of the 
statements sprinkled throughout the body of [2] that rely on Appendix B for justihcation 
are completely without merit. 


Appendix C: Appendix C purports to show that something is wrong with our second 
example given in [ 6 ]. (Our Erst example was criticized in subsection 4.2—see the 
discussion above.) For our second example, we simply wrote down a one-parameter 
family of metrics 5 'ab(A) that satisfied all of our assumptions except Einstein’s equation. 
We then declared it to solve Einstein’s equation with matter source TabW — Gab(A)/87r 
(“Synge’s method” 123 !). As one might expect, the resulting Tab{X) does not satisfy the 
weak energy condition. We then calculated the effective stress-energy of backreaction, 
tab, and showed that it failed to satisfy the conclusions of our two theorems (tracelessness 
and the weak energy condition). This example is useful because it shows that, for the 
validity of our theorems, we cannot dispense with the hypothesis of the weak energy 
condition on Tab (A). 

The authors of |2] correctly say that we obtain our matter stress-energy Tab{X) by 
Synge’s method, which they call “case (i).” They then say that another approach would 
be to obtain a different Tab(A) from a conformally invariant matter Lagrangian, which 
they call “case (ii).” Of course, case (ii) has nothing to do with our work. In case (ii), 
the matter stress-energy tensor would be conformally invariant, and their eq. (C.18) 


ofW call the “distributional derivative. ”] 

+ [NOTE ADDED: Despite the improved tone of Appendix B, the fundamental misunderstandings 
displayed in m remain present in m- For example, the first sentence of the paragraph following the 
paragraph containing eq. (B.20) in |I]/ asserts that “w-lim VaVf, 7 cd = 0 is not a consequence of the 
GW hypotheses (ii), (iii), (iv) but a strong a priori assumption ...” This assertion is false, since 
w-limVaV{, 7 cd = 0 is a trivial consequence of (A.l) and condition (ii).] 
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would hold. They then say, “Let us look more closely at Green and Wald’s choice (i): 
using (C.18) ...” They then proceed to derive contradictions. But eq. (C.18) applies 
only to case (ii), not to case (i). [NOTE ADDED: In |I]/, the reference to eq. (C.18) 
has been removed. Nevertheless, the association (made below eq. (C.22)) of the weak 
limit ofTab{g{l)) (which, obviously, is equal to Tab{g{l))) with Tab{0) makes no sense to 
us.] 

Since it is logically impossible for a metric and resulting stress-energy tensor 
obtained by Synge’s method to fail to be a solution to Einstein’s equation, it is logically 
impossible for our example to be wrong unless we made a computational error in 
calculating tab- We do not believe that we made any such error, and no evidence of 
such an error is presented in Appendix C. 
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